Hall Effect in spinor condensates 
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We consider a neutral spinor condensate moving in a periodic magnetic field. The spatially 
dependent magnetic field induces an effective spin dependent Lorentz force which in turn gives rise 
to a spin dependent Hall effect. Simulations of the Gross-Pitaevskii equation quantify the Hall 
effect. We discuss possible experimental realizations. 
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In 1879, Edwin Hall discovered that charges in a thin 
metallic plate exposed to a perpendicular magnetic field 
accumulate transverse to applied electric and magnetic 
fields. This Hall effect is used intensively to character- 
ize metals and semi-conductors. Modern interest in Hall 
effects started with the experimental discovery of the 
quantum Hall effect in two-dimensional electrons gases 
(2DEG) and the seemingly disconnected notion of Berry 
curvature [l[ and non abelian gauge field 0. Currently 
these notions explain the integer and fractional quantum 
Hall effect in 2DEG, the anomalous Hall effect (AHE) 
in ferromagnetic metals or semi-conductors and the spin 
analogue known as spin Hall effect 

Three mechanisms contribute to the AHE. Two of 
them, the side-jump [1] and the skew scattering [f| , result 
from carrier scattering off impurities while the Karplus- 
Luttinger term Q, also known as Berry curvature con- 
tribution, is related to spin-orbit coupling or non trivial 
magnetic order 0] in the band structure. The latter con- 
tribution has attracted considerable theoretical attention 
although it is very difficult to study experimentally in 
ferromagnetic compounds because of the additional con- 
tributions from impurities or defects. 

The realization of Bosc-Einstcin condensation in mag- 
netic or optical trap opens up the possibility of research 
at the border between atomic and condensed matter 
physics. Analogies between condensates in rotating mag- 
netic traps and electrons in strong magnetic fields al- 
low studies of the integer and fractional quantum Hall 
states |H . Recently, the realization of spin dependent 
optical lattices in combination with cold atomic gases (lol | 
offers the possibility to study spin-dependent transport 
phenomena. The neutrality of these gases, i.e. the ab- 
sence of a classical Lorentz force and the absence of im- 
purity scatterin g, a llows us to study the spin Hall ef- 



fect 



and the Berry curvature contri- 
I in a clean and controllable cn- 



bution of the AHE \U 
vironment. So far, the single particle approximation is 
often used to explain the AHE and SHE in cold atomic 
gases [III 12, EH- Although the single particle approx- 



scription of the Hall effect. 

Magnetic microtraps [l6| might also be used to study 
the AHE and SHE in Bose condensates. In this case, 
the Zeeman effect leads to a coupling between the spin of 
the atoms and the magnetic field of the microtrap. Var- 
ious geometries of magnetic field distributions including 
one-dimensional magnetic lattices have been explored ex- 
perimentally [IH , and two-dimensional magnetic lattices 
can be created using array of magnetic cylinders or 
magnetic dots for instance. 

In this letter, we investigate the Hall effect in spinor 
condensates where the magnetic field is created by a two- 
dimensional array of magnetic cylinders. The motion of 
the condensate through the magnetic lattice generates a 
non abelian gauge field that acts on the condensate as an 
effective spin dependent Lorentz force [l8|. Similarly to 
the case of fermions we first study the case where 
the two-body interactions are negligible. Then we study 
numerically the effect of the velocity of the condensate 
and of the non linear interactions on the Hall Effect. Fi- 
nally we propose an experimental setup where this effect 
can be observed. 

Let us consider a quasi-two-dimensional spinor conden- 
sate of spin F moving in a magnetic field B(r) created 
by a 2D array of magnetic cylinders. The Hamiltonian 
describing this system is (2(| 21, 22| 



U = j d 2 r^ a (r) 
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imation can be used to describe qualitatively the AHE 
in multicomponcnt condensates, two-body interactions 
should be included in the theory for a quantitative de- 



wherc ip a ( r ) with a = —F,...,F are the spin compo- 
nents of the quantum field operator ^(r), and V 2s is the 
projection operator that projects the spin state of the 
atoms pair into the total spin state 2s . F l=x,y > z are ma- 
trices of spin F and the constants a,2 S are the s-wave scat- 
tering lengths of the total spin channel 2s = 0, 2, . . . 2F. 
The confinement in the z direction appears in the Hamil- 
tonian via the oscillator length a z = y/h/Muj z , uj z = 
27r/ 2 where f z is the frequency of the trap along the z 
direction and M the mass of the atoms. The scalar poten- 
tial V(r) describes any external trap potential. Following 
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Ref. we apply a rotation T(r) of the quantization 
axis along the direction n(r) of the magnetic field B(r) 
to fTJ). The transformed Hamiltonian becomes 



|B(r)|K 

s=0 



2nh a,2s 
Ma z 



.A_(V r -iA s (r)) 2 + y(r)5 



(V 2s ) af}l s4(r)^(r) Mr), (2) 



where A 9 (r) = -iT+ (r)V r T(r) = A^(r)F* with i = 
x,y,z is the spin dependent gauge field resulting from 
the rotation of the spin axis. Explicit expressions for 
A*(r) (i = x,y,z) can be found in Ref. [23]. The two- 
body interactions are not affected by the gauge transfor- 
mation because they are short ranged and the projection 
operators are invariant against any rotation of the spin 
axis of both particles. Excluding the two-body interac- 
tions, Eq. © is a generalization of Eq. (2) of Ref. [l9[ for 
both fermions and bosons. The present work has two ad- 
ditional ingredients, the non linear interactions and the 
relaxation of the adiabaticity condition. The non abelian 
gauge field A s (r) contains terms that are proportional to 
F x and F y that induce transitions between the different 
spin-polarized states. 

Approximate analytical calculations: Let us initially 
assume that the spin-flip transitions can be neglected, 
i.e., that the spin adiabatically follows the magnetic field 
direction n(r). This approximation will be lifted later 
on. In contrast to F = 1/2, the expansion of Eq. @ 
gives rise to two spin-flip terms that are proportional to 
A g (r) and A s (r) • A s (r). The adiabatic approximation 

h ft 2 

is valid when Vi = < 1 and v 2 = 9rT ,^ A . < 1 



are fulfilled. The constant v is the group velocity of the 
particles, £ is a characteristic length of variations of n(r) 
and A z is the Zeeman splitting. The two-body inter- 
actions are included in Eq. [2] through a spin-spin inter- 
action term that can be decomposed into a scalar and 
a spin-dependent contribution [2jJ. The typical energy 
scale of the scalar inter-atomic interaction component is 
of the order of magnitude of the Zeeman splitting. The 
spin-dependent part of the two-body interactions, on the 
other hand, can be neglected because it is two to three 
orders of magnitude smaller than the scalar contribution. 
For this semi-classical analysis, we neglect the two-body 
interactions and assume that the spinor condensate is in 
a spin polarized state m E [— F, ...,F]. The effective 
Hamiltonian can be written as 



_h?_ 
2M 

m\B(v)\]ip m (r). 



d 3 r^ m (r) 



* A K r )) : 



V(r) 
(3) 



There is a difference between bosons and fermions. For 
fermions, all spin polarized states are affected by the 



gauge field A|(r) while for bosons, the spin polarized 
state m = is not affected by the magnetic texture. So 
for m ^ the field b a (r) = V x A*(r) = eijknid x njd y rik 
acts on the spin polarized state as an ordinary mag- 
netic field on a fictitious charge q = me (— e is the elec- 
tron charge), and then gives rise to an effective Lorentz 
force [IH that can induce a Hall effect. 

For a qualitative analysis, we consider the case where 
the average value of the field (b 3 (r)) s over one unit cell 
acts on the dynamics of the condensate. Firstly, the aver- 
age value (b s (r)) g = iirn/S with (n e Z) and S the sur- 
face of the unit cell, is quantized and depends only on few 
geometrical parameters such as the type and the period 
of the lattice. This value can be modified by applying 
a small constant magnetic field on top of the magnetic 
distribution Secondly, the spatial dependence of 

the field b s (r) and the potential |B(r)| are treated phe- 
nomenologically as causing momentum scattering that 
is characterized by the relaxation time r. We consider 
moreover that the condensate is under the influence of a 
force F = F x where F is the amplitude. Such forces 
can be induced by gravity for instance. The dynamics of 
the center of mass of the spin polarized state is described 
by: 



M- 



dt 



FqX — m \ 



h M x (b 



g/s 



(4) 



where (v) is the center-of-mass velocity of the polarized 
state. The solutions of Eq. (UJ) are simple and the wave 
packet trajectories go from closed orbits when there is no 
scattering and no force to open orbits when scattering 
and/or external force are present. One can also show 
that the quantity a xy = (v y ) / Fq is similar to the Drudc 
formula of the Hall conductivity component. The Hall 
effect is stronger when the spin F increases or the surface 
of the unit cell decreases. 

Numerical simulations: A quantitative study of the 
Hall effect in this systems can be done by solving Eq. (fTJ 
numerically using mean field theory. The order parame- 
ters of the spinor condensate, which are the expectation 
values of the field operators ip a {r), obey the time depen- 
dent Gross-Pitaevskii equation (TDGP) obtained from 
Eq. (TTj) by substituting i/j a {r) by its average. We solve the 
TDGP equation using an explicit Rungc-Kutta method 
of order six coupled to the Fast Fourier transform for the 
spatial derivatives. For the magnetic field distribution, 
we consider a square lattice of magnetic cylinders. 

The initial state of the system is obtained in the fol- 
lowing way. We first calculate the ground state of a 
spinor condensate of spin F = 1 composed of 2000 atoms 
of rubidium [24| in an harmonic trap with lu x = u> y = 
2-7T x 50 Hz (uj z — 2-7T x 1000 Hz). Then the condensate 
is shifted to one side of the trap and the wave-function 
is multiplied by a phase factor exp(ivx) where v is the 
group velocity. At t — the condensate is released from 
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FIG. 1: (a) Transverse position as a function of longitudinal 
coordinate of the center of mass of a F = 1 condensate for 
t^a, < 5 and different values of the initial group velocity v 
(in units of a x u> x ). (b) Time dependence of the ratio of aver- 
age transverse and longitudinal velocity. The initial state is 
prepared in the ferromagnetic state m = 1. 
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FIG. 2: Time dependence of the ratio (v v ) / (v x ) calculated 
for different initial polarizations and fixed initial velocity v = 
6. The vertical lines indicates the separation between the 
adiabatic regime where spin flip contributions are negligible 
from the regime where they dominate the dynamics of the 
condensate. 



the harmonic trap and is able to move freely through the 
magnetic lattice. 

To characterize the Hall effect, we calculate the ratio 
between the transverse and longitudinal velocity 



y 1 m 



/<«*> 



'.vy 



(5) 



where m indicates the initial polarization of the conden- 
sate at t = (in the global quantization axis which is 
along the z axis) and {vi =x . v ) is the average value of the 



velocity operator. a xy and a xx are the transverse and lon- 
gitudinal " conductivities" . This ratio corresponds phys- 
ically to the tangential of the angle between the electri- 
cal field and the velocity in the case of electronic sys- 
tems. For convenience, we express all quantities in re- 
duced units x — > a x x with a x — \J ^ ~ 1-59 yum and 
the energy e — » ehw x . Finally we choose the following 
set of parameters of the magnetic field. The period of 
B(r) is p = 7.5 a x , the amplitude is B = 10huj x and the 
distance separating the condensate from the surface of 
the magnetic array is h = 1.5 a x . 

The dynamical behavior of the TDGP equation is visu- 
alized in Fig. [T] which represents the ballistic trajectory 
of the center-of-mass of the Bose condensate (fig. [TJa) 
and the Hall angle (fig. [T]-b) for different values of the 
initial velocity. Fig. [l]-a shows that the trajectory of the 
condensate prepared in the state m = 1 is curved in the 
y direction because of the action of the gauge field on the 
dynamics of the condensate as qualitatively described by 
Eq. [3] Fig. [l]-b represents the time dependence of the 
ratio (v y ) m=1 / (vx) m= i- It shows pronounced oscilla- 
tions at early stages of the evolution that can be under- 
stood qualitatively by considering the sign of the effec- 
tive Lorentz force that is acting on the condensate. Since 
the size of the condensate is smaller than the size of the 
unit cell, part of the condensate will experience a positive 
Lorentz force while its complementary part will experi- 
ence a negative Lorentz force. Due to the distribution 
of particles in the condensate and the different topology 
of the regions where the field b a (r) is positive or nega- 
tive, the Lorentz force changes sign when the condensate 
is moving through the lattice and the y component of 
the group velocity oscillates with time. The oscillations 
of (v y ) =1 / (v x ) m=1 quickly disappear when t increases 
because (i) the condensate expands and (ii) the magneti- 
zation is not conserved so the mean field state describing 
the condensate tends to a state with cquipartition in pop- 
ulations (n±ifl = 1/3). Numerical simulations show that 
these oscillations remain when the size of the unit cell is 
smaller than the size of the condensate, as it should be 
expected. 

The dynamical properties of the system can be sep- 
arated into two different regimes; the adiabatic regime 
where the spin polarized state follows adiabatically B(r) 
as described by our approximate analytical calculations 
and a second regime that is dominated by the spin flip 
terms. The adiabatic regime can be identified at early 
time of the evolution of the condensate by calculating 
the ratio (v y ) / (v x ) for the spin polarized states m = 1 
and m = — 1. In this case, Eq. ([3]) shows that the two an- 
gles should differ by a sign which can be seen on the left 
part of Fig. [2] This behavior is not peculiar to Bose con- 
densates and should be observable with fcrmionic atomic 
gases or electrons gases provided that the gas is initially 
fully polarized and there is no disorder. As expected from 
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FIG. 3: (upper panel) schematic dependence of the Hall con- 
ductivity when a constant magnetic field is applied on the top 
of the magnetic lattice. The abrupt changes in (%)-, / {vx)± 
are related to the flux of b 9 (r) which changes sign at a pe- 
culiar value of Bo- (lower panel) numerical results for the 
ratio (v y } 1 / {v x ) 1 (5u^ 1 ) as a function of the magnetic field 
amplitude. 



amplitude of the constant magnetic field increases. 

To experimentally verify this theory, we propose to im- 
age the motion of a spin polarized condensate during its 
evolution in a microtrap with an additional periodic 2d 
magnetic field and deduce the ratio (v y ) / (v x ) by direct 
measurement of the velocity. The magnetic lattice can 
be created using arrays of magnetic cylinders or by pat- 
terning lattices of interconnected current loops. The con- 
densate can be accelerated through gravity by tilting the 
microtrap. 

To conclude, we have shown that spinor condensates 
in a magnetic lattice can be used to explore the geomet- 
rical phase contribution of the anomalous Hall effect. We 
found that this contribution is characterized by an adi- 
abatic regime which is valid in the short time limit and 
a long time regime where spin-flip processes dominate 
the condensate dynamics. We also demonstrated that 
the AHE is strongly affected by an additional constant 
external magnetic field. 

The authors thank M. Snoek, B. Canals, C. Lacroix, 
V. Dugaev, P. Bruno and J. Dalibard for discussion. This 
work was supported in part by Forschergruppe FOR 801 
of the Deutsche Forschungsgemeinschaft. 



Eq. ((3]), there should be no Hall effect for the polar state 
(m = 0) which can also be observed in Fig. [5J Fig. [2] 
shows that (v y ) / (v x ) Q of an initial polar state fluctu- 
ates around zero. These fluctuations are not affected by 
the spin-dependent part of the two-body interactions but 
are related to the spin-flip contributions of the gauge field 
A s (r). At longer time (right part Fig[5]), the dynamics 
of the condensate is governed by the spin-flip terms and 
the description of the Hall effect in terms of the adiabatic 
approximation in Q is not valid anymore. 

Finally, we study the influence of a small constant mag- 
netic field of amplitude Bq applied on top of the magnetic 
lattice in the direction perp endicular to the plane of mo- 
tion. Contrary to Ref. [19( which consider the electronic 
case, there is no classical Hall effect because the atoms 
are neutral. Therefore the features of the anomalous Hall 
effect are easily observed in Bose condensates. In the 
scmiclassical limit, i.e. when the spatial dependence of 
the gauge field is neglected, the variations of the Hall an- 
gle are abrupt because the average value of b a (r) can only 
change at peculiar values of Bq which depends on the lat- 
tice. Such variations are schematically represented on the 
upper panel of Fig|3] The general behavior described in 
the upper panel of Fig|3] can be reproduced numerically 
using the full time-dependent CP equation but the vari- 
ations of the ratio (v y ) / (v x ) are continuous rather than 
abrupt (lower panel of Fig|3]) . The numerical simulations 
also show that the Hall effect decrease rapidly when the 
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